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Abstract 

An experimental investigation of the effects of main- 
stream turbulence, mainstream swirl and non-symmetric mass 
addition has been conducted for the isothermal mixing of 
multiple jets injected into a confined rectangular crossflow. 
Jet penetration and mixing in the near field was studied using 
planar Mie scattering to measure time-averaged mixture frac- 
tion distributions. Orifice configurations were used that were 
optimized for mixing performance based on previous experi- 
mental and computational results for a homogeneous ap- 
proach flow. Mixing effectiveness, determined using a 
spatial unmixedness parameter based on the variance of the 
mean jet concentration distributions, was found to be mini- 
mally affected by free-stream turbulence but significantly 
influenced by the addition of swirl to the mainstream. The 
results for non-symmetric mass addition indicate that the 
concentration distribution of the flowfield can be tailored if 
desired. 

Nomenclature 

Aj orifice area 

A cross-sectional area of mainstream duct 

m 

at injection location 

C (S/H)VJ 

c av MR/(1 + MR) 

C d orifice discharge coefficient 

d orifice diameter 

DR P/Pn, 

H duct height = 5.08cm 

H eq effective duct height (see page 3) 

J jet-to-mainstream momentum-flux ratio 

= (pjVj 2 )/(p m V m 2 ) 

= (MR) 2 / [(DR) (Aj/AJ 2 (C d ) 2 ] 
mj mass flow of the jet 

m m mass flow of the mainstream 

MR nij/m m 

pj density of the jet 
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p m density of the mainstream 

S spacing between corresponding points of 

adjacent orifices, e.g. midpoints 
U s spatial unmixedness parameter 

Vj jet velocity = mj / (pj AjC d ) 

V m mainstream velocity = 304.8 cm/s 

V 1 root-mean-square of V m 

x downstream coordinate, x = 0 at the 

leading edge of the orifice 
y cross-stream coordinate (horizontal) 

z cross-stream coordinate (vertical) 


Introduction 

Crossflow mixing is convenient and efficient for a 
wide range of applications. Generally the objective is to 
rapidly obtain a homogeneous mixture of the injectant and 
mainstream. The appeal of the cross-stream configuration is 
its effectiveness and simplicity. The degree and rate of the 
mixing process is especially important in combustion applica- 
tions since exhaust composition directly depends on mass 
transfer and reaction kinetics. Whereas kinetics are difficult 
to control, the mixing process is easily affected by any number 
of parameters and optimization of that process for combustor 
design has been the topic of several recent investigations. 1 * 12 

In this paper the effect of the inlet flow condition of 
the mainstream on overall mixing performance in a rectangu- 
lar duct will be discussed. The downstream axial distance of 
interest is on the order of the duct height. Optimum configu- 
rations of jets, that is orifice size and spacing, have been 
previously identified using isothermal mixing experiments 5 12 
and CFD analysis 1 * 4 . Those results are based on a homoge- 
neous mainstream and jets that have initially symmetric 
velocity profiles injected perpendicular to the mainstream. 
However in practice the initial conditions of both flows are 
likely to be less controlled. Flow interactions caused by duct 
geometry and heat release, to mention a few sources, may 
affect the overall mixing process. The purpose of this inves- 
tigation is to compare the effects of non-ideal initial condi- 
tions on the development of jets in crossflow to previous 
results. 



Experimental 

Figure lisa schematic representation of the appara- 
tus which consists of 3 ducts of rectangular cross section. 
Sector width is 30.48cm. The inner duct height (H) is variable, 
but for these experiments was set at 5.08cm at the plane 
intersecting the row of orifices. The outer ducts which supply 
the jet fluid are H/2 in height. The inner and outer ducts are 
separated by removeable 3.18mm thick flat plates which 
contain the orifice configurations. For all tests near room 
temperature air was metered by 3 independently controlled 
critical flow orifices to each of the ducts. The baseline 
maximum variation in the mean approach velocity of the 
mainstream was experimentally determined to be 6% with a 
turbulence level of 1.3%. 

The effect of non-symmetric mass addition was 
studied by installing different orifice configurations on the top 
and bottom walls. The data were collected with the baseline 
mainstream flow condition. Subsequently perforated plates 
were placed in the mainstream to generate different levels of 
turbulence intensity (VW ). The plates were located 1H 
upstream of the plane of the orifice centers. Swirl was 
investigated by using a 4-swirler bulkhead placed 1.5H up- 
stream of the plane of the orifice centers. The flowfield 
induced by the perforated plates and by the swirlers was not 
experimentally characterized. 

Mie scattering was the primary diagnostic used to 
optically measure jet mixture fraction distributions in planes 
perpendicular to the duct axis (y-z plane). The planar digital 
imaging technique (for more detail see Ref. 13) relies on 
marking one of the flows with an oil aerosol (|im sized 
particles). A light sheet (0.51mm thick) created using a 2W 
argon-ion laser and a rotating mirror is used to illuminate a 
plane in the flowfield after passing through a window in the 
side wall of the test section. An image intensified thermo- 


electrically cooled CCD camera, located inside the duct 61cm 
downstream of the orifice center, was focused on the illumi- 
nated plane (end-on view). The camera was programmed to 
make exposures coincident with the sweep of the beam 
through the flow field. In these experiments the image was 
digitized in a 516 x 86 pixel format (pixel size = 0.58mm x 
0.58mm x 0.51mm) and sent to a computer for storage. The 
scattered light intensity is proportional to the number of 
particles in the measurement volume. If only one of two 
streams is marked, the light intensity of the undiluted marked 
fluid represents mole fraction unity. 

In this study, since the effect of different top and 
bottom jet mass flows was of interest, the mainstream was 
marked to avoid the more tedious setup of different marker 
concentrations for each of the jet flows in the non-symmetric 
mass addition study. The same setup was also used for the 
symmetric mass addition configurations. The data reduction 
procedure was used to reverse the marker location so that the 
intensity distributions presented are of the jet fluid. All of the 
reported measurements are mean concentration distributions 
for 15 second time-averaged measurements. 

Test Configurations 

The orifice configuration shown end-on in Fig. 2 was 
used to study the effect of unequal mass addition from oppos- 
ing sides of the duct, i.e. non-symmetric mass addition. Only 
the central 2H section of the 6H wide duct is shown. The top 
plate consists of circular H/d = 2.67 diameter orifices spaced 
S/H = 0.50 on center while the bottom plate consists of circular 
H/d = 4.00 diameter orifices also spaced at S/H = 0.50 on 
center. Both of these are optimum orifice configurations for 
a J of 25 as defined by the correlating expression for directly 
opposed inline orifices 2,6: 

c = (s/h)Vj (i) 



Figure 1 : Experimental Configuration used to Measure Planar Concentration Distributions 
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Figure 2: End-on view of the Non-symmetric Mass 
Addition Orifice Configuration. 


The trailing edge of each orifice was inline top-to- 
bottom as shown in Fig. 3. Mainstream flow is from left to 
right and the location of the downstream measurement planes 
are shown from the leading edge of the larger orifice, if the 
orifices on the top and bottom are different sizes. 



orifices on both side of the duct (test D and E) and then circular 
H/d = 4.00 diameter orifices on both sides of the duct (test F 
and G). In test E and F the opposing J values, and therefore 
mass-flow from each side, was held constant, while in test D 
and G the value of J was increased on the bottom row so that 
about twice the mass flow was introduced from the bottom vs. 
the top. 


The concept of a plane of symmtery between the top 
and bottom jets was utilized by Cox 14 and later Wittig et al. 15 
to allow empirical correlations to be developed for opposed jet 
injection using single sided injection data. Four different 
formulations that can be used to split the duct with an imagi- 
nary wall at the computed effective duct height (H^) are listed 
in Table 1 and defined below: 

The first formulation is a momentum-flux ratio balance: 
|H eq / H lj, l op = [ J W ([J] «op + [J1 bonom) 

The second formulation is based on Eq. 1, 

(where, C = (S/H)VX) : 

[H eq /H] Ctop =C lop /(C top + C bonom ) 

(note that (S/H), op = (S/H)^,,^ in all the experiments reported 
herein so this formulation cannot be separated from one that 
includes only VT ) 

The third formulation is a mass-flow ratio balance: 

[H«/H] M R.,op = MI V (MR-op + MR boaom > < 4 > 

(note that this gives the same answer as Wittig, et al. 15 when 

DR . 0 p = DR bo«o m and C d. .op = C d. bottom Since 

MR= (VT) (VDR)(C d ) (Aj/AJ) 

The fourth formulation is the one originally proposed by Cox 14 
based on a momentum balance: 

[H^flcox, top = KAj) (J)W (KAj) Wtop+KAjXJ) W < 5 > 


( 2 ) 

(3) 


Figure 3: Top view of the Non-symmetric Mass 
Addition Orifice Configuration. 


It follows that: 

[lyHUotn =l-[H eq /H] top 


( 6 ) 


Table 1 shows the test conditions for the 7 data sets 
used to illustrate the effect of non-symmetric mass addition. 
In the first three tests the non-symmetric orifice configuration 
shown in Fig. 2 is used. In test A, J is constant but the top vs. 
bottom mass- flow rate varies due to the different orifice areas. 
In test C, the top vs. bottom mass-flow rate is equal while J 
varies. In test B, both mass- flow rate and J vary. 

For comparison, the geometric configuration was 
then made symmetric by placing first circular H/d = 2.67 


In Table 2 the pattern (at half-scale) and the geomet- 
ric layout of the turbulence generating perforated plates used 
to study the effects of turbulence level on mixing performance 
are shown. Although experimental characterization was not 
performed, the level of turbulence intensity (V7V m ) indicated 
was estimated from the correlations found in Ref. 16. Four 
tests were performed with an orifice configuration consisting 
of circular H/d = 2.67 diameter orifices spaced S/H = 0.5 on 
center. The top and bottom plates were equil valent (Test E in 
Table 1). Data were only collected at x/H = 0.6875. 


3 




Table 1: Test Conditions used to Illustrate the Effect of Non-symmetric Mass Addition 
(note that H = 5.08cm in this investigation). 



Table 2: Turbulence Generating Perforated Plates 
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The effect of upstream swirl was investigated by 
using a 4-swirler bulkhead. A section of the bulkhead is 
shown in Fig. 4. The swirl is introduced through a 0.95cm 
annulus of inner radius 2.86cm. All of the mainstream flow 
passed through the annular passage and was swirled at 35 deg 
by the 1 8 vanes in each of the fixtures. Swirl is imparted in the 
counter-clockwise direction facing downstream. The bulk- 
head was located 1 .5H upstream of the jet injection location. 
The resulting flowfield was not experimentally characterized. 



Figure 4: End-on View of the Swirler Bulkhead. 


Analysis 


In a two-stream mixing problem the fully mixed 
concentration is defined by the jet-to-mainstream mass-flow 
ratio. A measure of the mixing rate can be obtained by 
comparing the jet mixture fraction distribution at any down- 
stream plane to the fully mixed value. In Ref. 1 1 the authors 
discussed a measure of unmixedness based on the variance of 
the concentration distribution, defined as spatial unmixedness: 


where, 



Cavg( 1 " Cavg ) 


(7) 


Ci 

Cavg 


c var ” ~r £ ( c i" c avg) 


= spatial concentration variance 
time-average concentration at a pixel 
fully mixed concentration 


U s = 0 corresponds to a perfectly mixed system, and 
U s = 1 a perfectly segregated system. The denominator is the 
maximum concentration fluctuation that can occur at the 
specified fully mixed concentration. Normalizing by this 
factor allows U s values to be compared regardless of the jet to 
mainstream mass-flow ratio of the system. Therefore, this 
parameter allows comparison of the relative mixing effective- 
ness of each configuration reported herein and comparison to 
other configurations with different mass-flow ratios. 



In many practical confined mixing applications the 
duct geometry is not symmetric about a plane that bisects 
opposing walls, for example in an annulus. However most 
experimental and computational studies have focused on 
configurations that are symmetric, i.e. opposing jets are equal 
in area and momentum-flux ratio (J) and therefore mass 
addition is symmetric. Correlations that predict optimum 
mixing performance based on these data sets, such as Eq.l, 
assume symmetry and generally utilize duct height (H) as a 
non-dimensionalization parameter. Besides being conve- 
nient, H is important since jet penetration is key to mixing 
performance. But what if mass addition is non-symmetric ? 
What is the correct value for the equivalent height ? 

The effect of non-symmetric mass addition is illus- 
trated in Fig. 5 by comparing the jet mass fraction distribution 
at four downstream locations for test A, B and C in Table 1 . A 
21-level color scale is used to represent contours ot jet mass 
fraction from 0 to 1 .0 (pure mainstream fluid colored red = 0 
and pure jet fluid colored dark blue = 1 .0). In each of the tests 
(A,B, and C) the top orifice configuration was identical and 
about twice the area of the bottom orifice configuration. In 
test A the opposing values of J are equal but the top vs. bottom 
the mass flow is not equal due to the different orifice area. The 
distribution shows a minimum on the duct centerline (H/2), 
which is not surprising since the jet trajectories were opti- 
mized for H/4 penetration using Eq. 1 . In test B the opposing 
mass flows are unequal and the values of J are unequal. The 
minimum in the concentration distribution is now nearer the 
top of the duct. In test C the opposing mass flows are equal but 
the values of J are unequal to obtain the mass balance through 
the different sized orifices. The values of J are even more 
mismatched than test B and the minimum in the concentration 
distribution is even nearer to the top than test B. 

The results of non-symmetric mass addition test A is 
compared to symmetric mass addition tests E and F in Fig. 6. 
All of these tests (A,E, and F) have the same value of J for the 
top and bottom orifice configurations. The minimum in the 
concentration distribution is seen to occur on the duct centerline 
in all cases. 

In Fig. 7 tests C, D, and G are shown. Each test was 
performed with J = 25 for the top orifices and J = 124 for the 
bottom orifices. Mass addition was symmetric in test C, but 
non-symmetric in D and G. However The minimum in the 
concentration distribution is seen to occur near the top of the 
duct in all cases. 
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Two of the formulations from Table 1 are plotted 
on Figs. 5 - 7 at the second downstream location. The solid 
line corresponds to the C formulation (Eq. 3) while the dashed 
line corresponds to the MR formulation (Eq. 4) (note that in 
Tests D,E, and F the formulations overlap). It appears that the 
better fit is obtained with the formulation based on Eq. 1, i.e. 
the value of H appears to depend on jet penetration which is 
determined by the product of S/H and the square root of J, 
rather than mass flow. Similar results were obtained by 
varying opposing values of J with inline circular H/d = 8 
diameter orifices spaced at S/H = 0.25 in Ref. 10 where both 
experiment and numerical results were reported. Although S/ 
H was not varied in either of these investigations the value of 
H «, appears to be independent of orifice diameter. 

Mainstream Turbulence 

The influence of several levels of freestream turbu- 
lence on the concentration distribution at x/H = 0.69 are 
compared in Fig. 8. The orifice configuration is symmetric 
with directly opposed circular H/d = 2.67 orifices spaced at 
S/H = 0.5. The value of J on both sides of the duct is 25 which 


is optimum for this configuration. The observed jet penetra- 
tions and mixing performance are similar even through the 
turbulence level is changed from 6 to 12 % using the perfo- 
rated plate geometries shown previously in Table 2. In fact the 
distributions are similar to the baseline trubulence level of 
1 .3% which was shown as test E in Fig. 6 and is shown again 
on Fig. 8. 

In order for mixing performance to be significantly 
affected, the penetration of the jet would have to be modified 
by the different turbulence levels. This is not evident in the 
data. It appears that the separation point and subsequent 
formation of the counter-rotating vortices is not affected by 
the turbulence level of the approach flow, although the data 
are only available at a single downstream position. This result 
is in agreement with the single jet experiments of Toften, et 
al. 17 where little change in vorticity was measured. It would 
appear that the vorticity generated by the velocity difference 
between the jet and mainstream dominates development of 
the flowfield, which is consistent with the results in Refs. 8 
and 9 where neither orifice shape nor downstream tabs had a 
significant effect on jet penetration or mixing. 


Test A 


Test B 


TestC 



Cavg = 0.55 Cavg = 0.58 Cavg = 0.62 

Figure 5: Effect of Non-symmetric Mass Addition 
(test conditions shown in Table 1) 


x/H = 0.375 
(trailing edge) 


x/H = 0.69 


x/H = 0.94 


x/H = 1.19 
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Mainstream Swirl 


The effect of the swirl on the concentration distribu- 
tion as a function of downstream position for circular H/d = 
3.70 orifices spaced at S/H = 0.36 at J = 50 is shown in Fig. 
9 (note that this is an optimum configuration). The effect of the 
approach flow is evident by the fact that immediately upon 
injection (trailing edge distribution) the top row of jets lean to 
the right and the bottom row of jets lean to the left. Farther 
downstream the "imprint" of the swirler is still evident in the 
distribution as periodic arrangement of jet fluid centered in the 
duct and corresponding to the location of the swirlers. The 
result suggests that the directional vector of the velocity is 
more important than the turbulence scale. Further work is 
required to characterize the approach condition and determine 
whether there is a significant effect on Eq. 1 . 

In Fig. 10 both the turbulent intensity and swirl 
spatial unmixedness curves are compared for the data in 
Figs. 8 and 9. Unfortunately the results of the turbulence 
level testing is only available at a single downstream 
position. The spatial unmixedness value for all of those 
points overlap and show up on Fig. 10 as a overlapping 
symbols at x/H = 0.688. The results of Test E are also 
plotted to compared optimized non-swirled performance to 
the swirled inlet condition. Although the actual test con- 
ditions are not identical, i.e. orifice configuration and J 
values, the performance should be similar since both orifice 
configurations are optimized using Eq. 1. It can be seen that 
the effect of the swirl is a slightly increased mixing rate 
whereas the turbulence level apparently has little effect. 






Conclusions 

• For non-symmetric mass addition an effective duct 

height (H ) based on Eq. 1 can be used to deter- 
mine the optimum mixing configuration for op- 
posed jet injection. 

• Mixing effectiveness for crossflow jets is minimally 

affected by free-stream turbulence. 

• Mainstream swirl can increase crossflow jet mixing 

rate. 
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